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Slow dynamics in random media: Crossover from glass to local- 
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Abstract. - We study slow dynamics of particles moving in a matrix of immobile obstacles 
using molecular dynamics simulations. The glass transition point decreases drastically as the 
obstacle density increases. At higher obstacle densities, the dynamics of mobile particles changes 
qualitatively from glass-like to a Lorentz-gas-like relaxation. This crossover is studied by density 
correlation functions, nonergodic parameters, mean square displacement, and nonlinear dynamic 
susceptibility. Our finding is qualitatively consistent with the results of recent numerical and 
theoretical studies on various spatially heterogeneous systems. Furthermore, we show that slow 
dynamics is surprisingly rich and sensitive to obstacle configurations. Especially, the reentrant 
transition is observed for a particular configuration, although its origin is not directly linked to 
the similar prediction based on mode-coupling theory. 
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Transport phenomena in inhomogeneous systems are of 
great importance in physics, chemistry, biology, and engi- 
neering [1] . Many biological systems and composite mate- 
rials consist of components of various sizes. The interplay 
between the broad range of length and time scales is essen- 
tial to their dynamical properties. Recently, slow dynam- 
ics and the glass transition of such systems has attracted 
much attention. One reason is that it is interesting to 
understand the effect of the generic disorder of inhomoge- 
neous systems on dynamic arrest [2] . Another reason is the 
extremely rich phenomenology that such systems exhibit 
near the glass transition point. For example, the colloidal 
suspensions with the short-range attractive potentials of- 
ten show a crossover from the glass transition at high den- 
sities to gelation at lower densities, which is triggered by 
the spatial disorder due to aggregation of the colloidal 
particles [3-6]. Other examples include peculiar glassy 
dynamics in binary colloidal mixtures with disparate size 
ratios [7-9], star polymer mixtures [10], or anomalous ion 
transport in silicate glasses [11-13]. However, the presence 
of multiple length/time scales have hampered elucidation 
of the origin of these interesting behaviors. For this reason, 
it is desired to study a simple model system where inherent 
complexities are pruned down as much as possible. The 
possibly simplest model is a mixture of mobile and immo- 



bile spherical particles. This system is a minimal model of 
spatially heterogeneous systems such as a fluid absorbed 
in porous media. It can also be seen as a model of a binary 
mixture where the characteristic time scales of constituent 
particles of each component are well separated. The slow 
dynamics of this model is interesting in its own right; in 
the dilute limit of immobile particles, dynamic arrest or 
the glass transition takes place at finite densities of mobile 
particles. The opposite limit where a single mobile par- 
ticle moves in a matrix of immobile particles is a classic 
Lorentz gas model. In this case, the mobile particle under- 
goes a localization transition at a finite immobile particle 
density, whose origin is purely geometric [14,15]. Recently, 
Krakoviack studied this model theoretically for arbitrary 
densities of mobile/immobile components by combining 
the replica method with mode-coupling theory (replica 
MCT or RMCT) [16-18]. Two important predictions have 
been made in his studies. One is the crossover of dynam- 
ics from the glass to localization transition. When the 
volume fraction of mobile particles, <f> m , is large and the 
volume fraction of immobile particles, <pi, is small, the 
system undergoes a glass transition, where the onset of 
slow dynamics is signaled by the discontinuous emergence 
of two-step relaxation, over all wavelengths, in the den- 
sity correlation functions. This transition is referred to as 
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Fig. 1: Dynamic phase diagram of a bidisperse soft-sphere sys- 
tem generated by the EM protocol. The dotted line represents 
the Liquid- Arrested (LA), line which is defined as the points 
beyond which the relaxation time r a exceeds 10 3 . 



Type B transition in the glass transition community [19]. 
As fa increases, the glass transition point fa\ decreases 
drastically. At even larger fa, dynamics near the transi- 
tion point qualitatively changes; one-step slow relaxation 
sets in at large wavelengths, where the amplitude of the re- 
laxation curve grows continuously and progressively prop- 
agates toward the shorter wavelengths as fan increases. 
This dynamics is called Type A. The crossover from Type 
B to A is especially interesting since it is ubiquitous in 
many heterogeneous systems such as colloidal gels [4, 6] 
and systems with disparate size ratios [9]. Another pre- 
diction of RMCT is the existence of the reentrant pocket 
in the low fan and high fa region, in which the dynam- 
ics of mobile particles is accelerated as fa increases. This 
counter-intuitive result was rationalized as being caused 
by dephasing of the dynamics of the caged particles due 
to multiple collisions with other mobile particles [16-18]. 

In this Letter, we investigate slow dynamics of mix- 
tures of mobile/immobile particles by extensive numeri- 
cal simulations. Our goals are twofold; first, we study 
the crossover from the glass to localization transition, ex- 
ploring the whole range of parameter space of {fa, fan)- 
Slow dynamics of mobile particles near the transition line 
is analyzed using the density correlation functions, non- 
ergodic parameters, and mean square displacement. We 
also calculate the dynamic susceptibility in order to probe 
the dynamic heterogeneities around this crossover. The 
second goal is to study the sensitivity of the dynamics of 
mobile particles to a geometry of random configurations 
of the immobile particles. To see this, we investigate two 
types of protocols to generate random configurations. In 
the first protocols which we refer to as Quenched- Annealed 
(QA), we prepare the immobile particles, let them run 
for some time steps and, after they are equilibrated, we 
quench their motions. Then, the mobile particles are in- 
serted into the void spaces of the matrix of immobile par- 



ticles. This QA protocol is appropriate to study dynamics 
in a porous medium. This is the protocol that has been 
adopted in RMCT analysis [16-18]. In the second proto- 
col which we call Equilibrated Mixture (EM), we first run 
all of particles until they are equilibrated and then freeze 
the motion of a fraction of the particles while keeping the 
others mobile. We design this protocol in order to mimic 
binary mixtures with disparate time scales [7-9, 20] and 
to study the dynamics of the faster (lighter) component, 
which would be analogous to that of mobile particles in an 
EM system. In our study, we observe no sign of the reen- 
trant transition for QA systems, which is contrary to the 
RMCT prediction [16-18]. Surprisingly, however, we do 
find the reentrant transition for EM systems. The phys- 
ical origin of this reentrance is purely geometric and no 
direct connection with the RMCT prediction is inferred. 

We carried out molecular dynamics (MD) simulations 
for two systems. The first one consists of an equal num- 
ber of two types of particles. The reason to use the 
bidisperse system was to avoid crystallization in large 
fan and small fa regions. Our system is composed of 
N = Ni + N2 = 1000 particles in a cubic box with 
dimension L = V 1 ^ = 10.8 under periodic bound- 
ary conditions (PBC). They interact via the soft poten- 
tial v ab (r) = e(er Qb /r) 12 , where a ab = (a a + er 6 )/2 and 
a, b 6 {1,2}. The size and mass ratio were o^/fi = 1-2 
and milrtix = 2, respectively. The units of length, time, 
and temperature were taken as 01, yj rniaf/e, and e/ks, 
respectively. Ni (N m ) immobile (mobile) particles were 
chosen from the N\ + N2 particles for each simulation run. 
We used the volume fractions defined by fa = nNia^/QV 
and fa n = nN m a^ s /6V as the system parameters, where 
CTcff = ? 1 - 1 / 12 (^ Q b cr3 f) /4) 1 / 3 is the effective diameter of 
particles [21]. The states investigated here were a^ s = 0.5, 
0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.15, 1.2, 1.3, 1.4, and 1.45. The 
corresponding temperatures were T = 21.61, 10.42, 5.624, 
3.297, 2.058, 1.350, 0.992, 0.772, 0.651, 0.473, 0.352, and 
0.306, respectively. We controlled fa and fan by changing 
Ni, N m , and <7 c ff. The minimum number of mobile par- 
ticles was chosen as N m = 10. The time steps used were 
0.001 ~ 0.005. The second system we investigated is the 
monatomic hard spheres. This system was employed in 
order to clarify the origin of the configuration dependence 
of the dynamics without being obscured by the bidisper- 
sity and softness of the continuous potential. The system 
includes N = 1000 hard spheres with mass m and diam- 
eter a in a cubic box of volume V under PBC. a and 
\Jraa 1 /ksT were used as the units of length and time, 
respectively. The temperature was fixed as fcgT = 1. The 
volume fractions fa and fan were controlled by changing 
Ni, N m , and V. For both systems, wc carefully checked 
the sample dependence of the observables throughout the 
study. 

We first performed simulations of EM systems of the 
bidisperse soft-sphere system to observe the dynamics for 
the whole (fa, m )-space. In Fig. [H the dynamic phase di- 
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Fig. 2: (a)-(c): The self part of the intermediate scattering 
function F a (k,t) at k = 2n. (a)-(c) corresponds to the cuts 
indicated in Fig. [T] From left to right, the densities (4>i,4> m ) 
increase along the arrows denoted in Fig. [1] Inset of (c) : k- 
dependence of the nonergodic parameter, (d): Mean square 
displacement for the cut (c) in Fig. [T] From top to bottom, 
the densities (fa,4>m) increase. 

agram is plotted as a function of fa and <j) m . This is drawn 
by calculating the self part of the intermediate scattering 
function (ISF) for mobile particles defined by F s (k,t) = 
N ™ ( E 4 ="i e l ^W-^°»). Here k is the wavevector, 
k = \k\, and fi(t) is the position of the i-th particle. We 
define the Liquid- Arrest (LA) line that delineates the liq- 
uid and arrested (glass or localization) phase in Fig. Q] 
as the points at which the structural relaxation time r Q 
reaches 10 3 . Here r Q is defined by F s (k = 27t,t q ) = 0.1. 
It is seen that, in the small fa regime, the LA line (or 
the glass transition points) drastically decreases as fa in- 
creases. This behavior is in accordance with the results of 
simulations [22-25] and RMCT [16-18]. This trend is sus- 
tained up to the small <f> m regime, beyond which a small 
reentrant pocket is observed. We shall discuss this reen- 
trance later in detail. 

In Fig. H] (a)-(c), we show F s (k,t) for the mobile parti- 
cles along the three cuts across the LA line, as marked by 
arrows (a)-(c) in Fig. [TJ We also calculated the collective 
part of the ISF, F(k,t) and found its behavior is qualita- 
tively the same as F s (k, t) except at small wavevectors, the 
details of which will be explained elsewhere [26]. Figure [2] 
clearly demonstrates the existence of two types of distinct 
dynamics depending on the immobile particle densities. 
As shown in Fig. Ufa), at small fa, F s (k, t) exhibits a typi- 
cal glassy behavior; the two-step relaxation with a plateau. 
This plateau or the nonergodic parameter (NEP) appears 
dis continuously over the whole wavelength range at an 



onset density, followed by a mild increase of the plateau 
heights as the densities (fa, cj> m ) increase along the arrow. 
This is typical Type B behavior. At large fa, on the othcr 
hand, the ISFs show a single step relaxation with a long 
tail. Its amplitude or NEP continuously grows as (fa, <f> m ) 
increases across the LA line, as seen in Fig. |2](c) . The k- 
dependence of the NEP is plotted in the inset of Fig.[2jc). 
NEP (calculated from F(k,t) / F(k,0)) emerges from zero 
continuously at small k, as (fa, 4> m ) approach the LA line. 
As (fa,4> m ) increase further, NEP grows and this trend 
propagates progressively toward larger k. These behaviors 
are the hallmark of Type A or the localization transition 
predicted by RMCT and also demonstrated by simula- 
tions for various spatially heterogeneous systems such as 
colloidal gels [6,9]. The behavior of the mean square dis- 
placement (MSD) defined by 8r 2 (t) = A^ x (ES Aff (i)), 
where Afi(t) = fi(t) — ^i(O), also qualitatively changes as 
the dynamics changes from Type B to A. It is well known 
that, in the Type B regime, MSD shows a plateau caused 
by the cage effect (the results not shown). On the other 
hand, at the low 4> m limit, where the system becomes a 
Lorentz gas, the anomalous subdiffusive behavior is ob- 
served at long times, which can be explained in terms of 
percolation theory [15] . The latter behavior is observed for 
the MSD in the Type A regime, as shown in Fig.[D(d). The 
MSD shows regular diffusive behavior in the liquid region, 
but near the LA line around (fa,fa n ) ~ (0.377,0.042), it 
becomes subdiffusive; Sr 2 (t) « t - 3 , the exponent of which 
is consistent with 0.32 predicted for a Lorentz gas [27,28]. 

We next investigated the nonlinear dynamic suscepti- 
bility or 4-point correlation function, X4(fc,t), for mo- 
bile particles [29-31], which is a good measure of the 
dynamic heterogeneities. X4(fc,t) is defined as the vari- 
ance of the fluctuations of the self part of the ISF by 
Xi (k,t) = N rn [(F s (k,t) 2 ) - (F s (k,t)) 2 }. Here F s (k,t) = 
(F s (k,t)). In the literature [30], N~ l cos(k ■ An(t)) 
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is used as the definition of F s (k,t). With this defini- 
tion, however, Xi(k,t) decays to 1/2 at t — > oo. On 
the other hand, as we demonstrate here, the peak of 
the dynamic susceptibility for the Type A regime grows 
more mildly than for the Type B regime. In order to 
demonstrate this suppression of the peak and thus the 
dynamic heterogeneities in the Type A regime without 
being obscured by a constant plateau at large t, we 
have adopted an alternative definition of F s (k, t) which is 
Nm 1 T,^sm{k\Af i (t)\)/k\Afi(t)\. Of course, both def- 
initions of F s (k,i) lead to the identical averaged value 
F s (k,t) = (F s (k,t)). With this new definition, x4M) 
decays to at t — ► oo. In Fig. [31 the time evolution of 
Xi{k,t) is plotted for the three states (a)-(c) in Fig.[TJ In 
the glass (Type B) regime at (a), X4(fc,i) exhibits behav- 
ior well-known for the bulk glass, i.e., a pronounced peak 
at the cv-relaxation time, whose height grows rapidly as 
the density increases, preceded by algebraic growth in the 
/3-relaxation regime. In the localization (Type A) regime 
at (c), on the other hand, X4(^j^) does not grow nor show 
a strong peak, even at long times. This implies that dy- 
namic heterogeneities play a minor role in the slow dy- 
namics in this regime. A similar behavior of X4(fc,t) has 
been reported for colloidal gels whose slow dynamics is 
caused by geometrical frustrations [32,33]. 

Finally, we examine the sensitivity of the above results 
to the configurations of immobile particles by comparing 
results for the Equilibrated Mixture (EM) with those us- 
ing the Quenched-Annealed (QA) protocol. Hereafter, we 
mainly focus on the moderately small <j) m regime, because 
little protocol dependence is observed at large <p m . Here 
we switch the system to the monatomic hard sphere sys- 
tem to keep physical interpretation from being obscured 
by the bidispersity and continuous potential of the pre- 
vious system. We checked that the results remain un- 
changed for the bidisperse soft-sphere systems. Dynamic 
phase diagrams of both QA and EM systems are plotted 
in Fig. |U Here the LA line is determined as the points be- 
yond which the MSD does not exceed 10 2 in a time within 
the simulation time t = 10 4 . First, from the simulation 
for the lowest <p m of our study (4> m = 0.026), we estimate 
the percolation point to be 4% « 0.25 for both QA and 
EM systems, which is close to 0.24 determined by Voronoi 
tessellation [24]. In the proximity of the MSD again 
shows subdiffusive behavior 8r 2 (t) k t 3 [27,28]. Details 
of these results will be explained elsewhere [26]. We here 
note that to check the finite size effects we performed the 
same simulations for N = 10000 and confirmed that over- 
all behavior is almost the same in both the QA and EM 
systems. 

As seen in Fig.^a), no reentrant transition is observed 
for QA systems, contrary to the RMCT prediction. The 
LA line monotonically decreases as 4>i increases, which is 
compatible with the recent numerical simulations for the 
related QA system [34]. This result is hardly surprising, 
considering that the slow-down of mobile particle dynam- 
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Fig. 4: Dynamic phase diagram of hard sphere systems with a 
random matrix generated with two different protocols; (a) QA 
and (b) EM protocols. 



ics is mainly due to the geometrical confinement by the 
immobile particles. Krakoviack has conjectured that the 
reentrant pocket predicted by RMCT is caused by the de- 
phasing of the trajectories of the caged mobile particles 
by collision with other mobile particles [16-18]. This ef- 
fect may be too small to manifest itself as reentrance, if 
any. On the other hand, surprisingly, the EM system has 
the reentrant pocket, as clearly seen in FigHJb) and also 
in Fig.l. In other words, the dynamics of the mobile par- 
ticles is accelerated as the number of mobile particles in- 
creases. The physical origin of the reentrant pocket for EM 
systems is clear. The random configuration of the immo- 
bile particles is generated by quenching their motion after 
both mobile/immobile particles arc equilibrated together. 
Therefore immobile particles adjust themselves so that the 
free volume of both components is entropically maximized 
and leaves more available geometrical pathways for mobile 
particles, which pushes the percolation point 4\ to larger 
values and thus results in faster dynamics of mobile par- 
ticles. Sensitivity of the percolation point on the proto- 
cols used to generate the configuration has been studied 
in several contexts [23-25]. However, its interplay with 
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the dynamics of mobile particles, especially the reentrant 
transition, has not been explored, to the best of our knowl- 
edge. It is worth mentioning that an analogous reentrant 
behavior has been found for a distinct system, i.e., a bi- 
nary mixture consisting of large/small particles studied by 
MD simulations [20]. It has been observed that diffusion 
of small particles becomes slower when the interactions 
between small particles are switched off [20]. We believe 
that this counter-intuitive observation is closely related 
to our finding of the reentrant pocket for EM systems this 
switch-off leaves small particles "invisible" from each other 
and the large ones feel as if there are fewer small particles 
around them. Therefore, the equilibrium configurations 
of the large (thus less mobile) particles would leave less 
void spaces and pathways for small particles than before 
the switch-off. This consideration and the results of the 
present study naturally lead to a speculation that separa- 
tion of time scales, rather than disparity of particle sizes, 
is a more essential element behind anomalous dynamics 
(such as reentrance) in heterogeneous systems. 

In summary, we have carried out MD simulations to 
study slow dynamics in random media for a broad range 
of parameter space. Two types of dynamics, A and B, 
have been observed and quantified in terms of various dy- 
namic quantities, which supports RMCT predictions qual- 
itatively. The dynamic susceptibility clearly demonstrates 
the minor roles of dynamic heterogeneities in the Type A 
regime. We did not observe any evidence of the reentrant 
transition for the system for which RMCT has been origi- 
nally applied. Contrary, we do observe the reentrant tran- 
sition for the system with configurations generated with 
a different protocol. However, the physical origin of this 
reentrance is not directly linked to the similar prediction 
based on RMCT. Our study suggests that anomalous dy- 
namics common in the heterogeneous systems might be, 
more often than not, explained well by a maximally sim- 
plified model such as the one in the present study. 
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